Compositional Inversion Symmetry Breaking in Ferroelectric Perovskites 
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Ternary cubic perovskite compounds of the form (A 1 /3A' 1 ^ 3 A"y 3 )B03 and A(B 1 /3B' 1 y 3 B' 1 '^ 3 )03, in 
which the differentiated cations form an alternating series of monolayers, are studied using first- 
principles methods. Such compounds are representative of a possible new class of materials in 
which ferroelectricity is perturbed by compositional breaking of inversion symmetry. For isovalent 
substitution on either sublattice, the ferroelectric double-well potential is found to persist, but 
becomes sufficiently asymmetric that minority domains may no longer survive. The strength of the 
symmetry breaking is enormously stronger for heterovalent substitution, so that the double- well 
behavior is completely destroyed. Possible means of tuning between these behaviors may allow for 
the optimization of resulting materials properties. 

PACS numbers: 77.80.-e, 77.84.Dy, 61.50.Ah, 81.05.Zx 
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In the last decade, the extensive use of first-principles 
theoretical methods to study ferroelectric perovskite ox- 
ides has led to a greatly expanded understanding of the 
ferroelectric and piezoelectric properties of this impor- 
tant class of materials. Theoretical investigations of elec- 
tronic, dynamical, and structural properties have been 
shown to be in good accord with experimental observa- 
tions for the simple ABO3 perovskites Q. Moreover, 
these studies provide microscopic insight into the fer- 
roelectric instabilities, their relation to the long-range 
Coulomb interactions g] , and related questions about the 
origins of the piezoelectric response H,Q] . 

However, the materials of most interest for techno- 
logical applications are generally not the simple ABO3 
perovskites, but solid solutions with stoichiometric sub- 
stitutions of A or B metal atoms. Examples include 
PZT (PbZi'i-^Ti^Oa), currently one of the most widely 
used ferroelectrics, and PZN (PbZni/3Nb2/303) and 
PMN (PbMg 1/3 Nb2/303) and their solid solutions with 
PbTiC>3, which have recently been shown to have enor- 
mous piezoelectric response in single-crystal form H . In 
fact, this class of materials shows great promise for the 
development of new materials having improved dielectric 
and electromechanical properties. Not only is there an 
enormous space of chemical compositions to explore, but 
it may also be possible to optimize the desired material 
properties by tuning the degree and type of compositional 
order for the desired application 

In particular, Eckstein || has recently suggested that 
the artificial atomic-layer growth of compositionally or- 
dered structures that break inversion symmetry might be 
especially exciting and fruitful in this regard. The result- 
ing asymmetry of the ferroelectric double-well potential 
in such a material suggests the prospect of qualitatively 
new behavior, e.g., "self-poling" materials with tailored 
piezoelectric or dielectric properties. 

In this Letter, we explore this exciting prospect by 
carrying out ab-initio theoretical calculations for several 



prospective model structures of this type. Specifically, we 
envision the artificial growth of materials of overall com- 
position (A 1/3 A' 1/3 A' 1 ' /3 )B03 or A(B 1/3 B , 1/3 B' 1 ' /3 )03, in 
which the three different cations alternate layer-by-layer 
along the ferroelectric direction as illustrated in Figs. |](a- 
b) for A-site and B-site modulation respectively. We will 
assume that such atomic-layer control will become possi- 
ble J§|, and investigate the energy landscapes and ferro- 
electric properties of the resulting materials. As will be 
shown below, we find that the asymmetry may easily be 
strong enough for self-poling to occur, even in the case 
of isovalent substitution. Moreover, we find a surpris- 
ingly strong variation of the strength of the symmetry 
breaking with the strength of the compositional pertur- 
bation, allowing for a very wide tunability of ferroelec- 
tric properties. Experimentally, such tunability might be 
exploited by growing alloy structures with a cyclic mod- 
ulation of the concentration variable. We thus find very 
strong motivations for the development of such materi- 
als, and hope that our work will encourage experimental 
efforts directed toward their synthesis. 

We adopt bulk BaTiC>3 as a prototypical parent com- 
pound and construct a series of model systems that al- 
low us to test for the effects of A-site vs. B-site and 
isovalent vs. heterovalent substitution. Specifically, we 
consider (Ba 1 / 3 Sr 1 / 3 Cai/3)Ti03 as an example of an 
"A-iso" system (isovalent substitution on the A site), 
Ba(Sci/3Ti 1 / 3 Nbi/ 3 )03 as an example of a "B-hetero" 
system (heterovalent substitution on the B site), and 
Ba(Ti 1 / 3 Zr 1 / 3 Hf 1 / 3 )0 3 as an intermediate "B-iso" case. 
We construct 15-atom supercells by tripling the primi- 
tive unit cell along the z direction and cyclically alter- 
nating the identity of the A or B atom, as shown in 
Figs, l(a-b). We assume that any ferroelectric order de- 
velops only along the z-direction, so that the material 
remains tetragonal, and only z displacements need be 
considered. (It may be possible to realize this situation 
by appropriate choice of lattice- mismatched substrate for 
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FIG. 1. Side views of (a) (Ai/sAi/gA'/^BOs, and (b) 
A(B 1/ /3B' 1 / 3 B / 1 ' / / 3 )03 structures. Energy vs. displacement 
along (c) the cubic force direction, (d) the line r = rx + wo£o, 
and (e) the line connecting the minima A and B. Solid 
(dashed) curves are for A-iso (B-hetero) systems. 



epitaxial growth, but it is not our purpose here to inves- 
tigate this possibility.) We also assume perfect control of 
layer-by-layer composition, resulting in the ideal stacking 
of Figs, l(a-b), and carry out our theoretical studies only 
at zero temperature. Much future work clearly remains 
to be done in relaxing these assumptions, but the study 
of simple prototypical systems is a natural first step. 

The ab-initio calculations are carried out using the 
Vanderbilt ultra-soft pseudopotential scheme ||] in the 
local-density approximation. Details of the pseudopo- 
tentials can be found in Ref. pX] . Good fc-point conver- 
gence is obtained using a (6,6j2) Monkhorst-Pack mesh, 
corresponding to the bulk (6,6,6) mesh [jfl|, and a 25-Ry 
plane- wave cutoff is used throughout. Because experi- 
mental lattice constants a and c are not available for the 
ordered compounds of interest, the calculated theoretical 
equilibrium values were usually used. However, in some 
cases the c/a ratio was fixed and only the cell volume 
was optimized, as detailed below. 

In the presence of the broken inversion symmetry, it 
is not always easy to locate both local energy minima 
corresponding to the two ferroelectric ground states. We 
have found that the following procedure is quite reliable 
for finding both minima, if they exist. First, we place 
the atoms at the ideal cubic coordinates and calculate 
the pattern of forces, thus defining a "cubic force direc- 
tion" £ c f in the 15-dimensional configuration space. We 
then use £ c f as the search direction for a line minimiza- 
tion, and denote the minimum along this line r = u c f £ c f 
to be X, as shown in Fig. 0(c). Next, starting from 
X, we carry out a second line minimization along the 
line r = r x + uq £o , where £o corresponds to the ferro- 
electric mode unit vector of bulk BaTiC>3. As shown in 




FIG. 2. Energy vs. displacement along the line 
r = rx + uq cjo in the A-iso system for several cell volumes. 



Fig. |l](d), two scenarios can be identified, (i) If two min- 
ima A' and B 1 are found along this line, then we carry 
out a steepest-descent minimization from each, obtain- 
ing two distinct local minima A and B. Fig. |l|(e) shows 
the perturbed double-well potential plotted along a line 
connecting points A and B (direction £) in such a case, 
(ii) If only one minimum is found, as for the dashed line 
of Fig. |l](d), then the double- well potential has been de- 
stroyed, and a subsequent steepest-descent minimization 
identifies the unique minimum. 

We first consider the application of this procedure to 
the A-iso system (Bai/3Sr 1 / 3 Cai/3)Ti03. The relaxation 
from the ideal cubic structure to the point X illustrated 
in Fig. 0(c) is found to be dominated by ionic size ef- 
fects: one sees a simple shift of the atoms toward smaller 
A-site cations and away from larger ones. Figure ||(d) 
then illustrates the minimization along r = r x + uo £o f° r 
one particular fixed cell volume (1343 a. u. 3 ), indicating a 
modest asymmetry between minima A' and B' . However, 
since it is well established that the presence or absence 
of a ferroelectric instability depends strongly on cell vol- 
ume in perovskites [ |To| , we have repeated this step for a 
series of fixed atomic volumes in Fig. (The curves in 
Fig. H could thus be roughly interpreted as correspond- 
ing to a series of related compounds having differing ten- 
dencies towards ferroelectric instability.) The results il- 
lustrate the typical evolution of the double-well poten- 
tial with increased ferroelectric tendency, from a slightly 
asymmetric single-well minimum (V=1123 a.u. 3 ) to a 
slightly asymmetric double- well potential (F=1463 a.u. 3 ) 
p"2[ . It happens that the theoretical equilibrium volume 
l^=1193a.u. 3 is a nearly borderline case; we find only a 
single minimum, but there are still two inflection points. 
In such a case, we modify the procedure of the previ- 
ous paragraph by carrying out a steepest-descent mini- 
mization from each inflection point, and by doing so we 
succeed in finding the two distinct minima and in 
Fig. |l|(e). For this case (V^=1193a.u. 3 ), we find a soft- 
mode amplitude of 0.3 a.u. and an average well depth (rel- 
ative to the saddle point S) of 0.48 mHa per 5-atom cell, 
fairly close to the corresponding values of 0.25 a.u. and 
0.43 mHa ju^, respectively, for tetragonal bulk BaTiC>3. 
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The energy difference between the twin weifs is modest, 
- 15 - 20% of the well depth. 

Applying the same approach to the B-iso system 
Ba(Ti 1 / 3 Zr 1 / 3 Hf 1 / 3 )03, we observe very similar behav- 
ior: a strong volume-dependence of the ferroelectric ten- 
dency, and a modest inversion asymmetry for both single- 
well and double-well volumes. In this case, the structure 
does not develop a double-well potential until the vol- 
ume is increased to 1468 a. u. 3 , but the degree of asym- 
metry is similar as for the A-iso case. Thus, we conclude 
that the choice of A vs. B site for an isovalent chemical 
substitution does not strongly affect the strength of the 
asymmetry or the qualitative behavior of the system. 

While the asymmetries may appear small in the case 
of isovalent substitution, they are large by one impor- 
tant measure. We define an "effective electric field" 
£ c g = AE/AP S , where AE and AP S denote the energy 
and polarization differences between the two local min- 
ima. If this quantity is larger than the coercive field of 
the material, it means that the thermodynamic prefer- 
ence for the energetically preferred minimum is strong 
enough to overcome the pinning of the domain walls and 
to spontaneously switch the material into a single-domain 
state. Using the Berry-phase approach jl3| to calculate 
AP S , we find £ G ff=90kV/cm at the equilibrium volume 
in the A-iso system. Since the typical coercive field of 
most perovskite ferroelectrics is closer to 15kV/cm, we 
thus arrive at the important conclusion that the symme- 
try breaking may easily be strong enough to cause the 
material to self-pole, even in the case of isovalent substi- 
tution. 

Turning now to the case of heterovalent chemical 
substitution, we find very different behavior in this 
case. Specifically, we consider the B-hetero system 
Ba(Sc 1 / 3 Ti 1 / 3 Nb 1 / 3 )03 in which the valence charges are 
+3, +4, and +5 on Sc, Ti, and Nb, respectively. Struc- 
tural optimization results in a lattice constant of a = 
7.60 a. u. and c/a = 3.023. The pattern of relaxation £ c f 
leading to configuration X suggests that the electrostatic 
interaction is a dominant effect. Specifically, we observe 
bucklings of the AO and BO2 planes that are consistent 
with a picture of static electric fields arising from the 
different B-atom valence charges. 

Most importantly, for the heterovalent case we find 
only a single minimum when searching from X along 
£0, as shown by the dashed curve in Fig. 0(d). Vari- 
ous alternate search strategies failed to identify a second 
minimum; all trial structures were found to relax back 
to a unique structural ground state. The absence of the 
second minimum was not found to be sensitive to the 
cell volume, as was the case for isovalent substitutions. 
Thus, we find that the symmetry breaking is enormously 
stronger than for the isovalent case, and it is clear that 
resulting behavior is of a qualitatively different type. 

In order to gain a better understanding of this be- 
havior, and in particular to track the disappearance of 
the secondary minimum, we have developed a model sys- 
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FIG. 3. Energy vs. displacement along the line 
r = vx + u £0 in the B-hetero system with continuous 6. 



tern to study the effects of "turning on" the heterova- 
lent symmetry-breaking perturbation gradually. We in- 
troduce a continuous variable S, and construct artificial 
atoms with fractional nuclear charges that deviate by 
±<5 from those of Ti (Z—22). Constructing a crystal 
out of a cyclic series of Ti— S, Ti, and Ti+<5 atomic lay- 
ers, we can continuously tune the system from ferroelec- 
tric BaTi0 3 (<5=0) to a full-fledged heterovalent system 
Ba(Sc 1 / 3 Ti 1 / 3 V 1 / 3 )0 3 (S=l). The energies calculated 

along the £0 direction from X are plotted for a series of 8 
values in Fig. |. As 5 increases from to 0.3, the asym- 
metry of the double well increases and the well depth 
decreases. At 5 — 0.4, the curve exhibits only a single 
energy minimum, signaling the transition from double- 
well to single- well behavior in response to the increasingly 
strong symmetry-breaking perturbation. To confirm the 
disappearance of the secondary minimum more directly, 
we also tracked its evolution as S was "turned on" in a 
sequence of small steps, using the relaxed structure at 
one S as a starting guess for the next. At 6 — 0.34, the 
minimum was confirmed to disappear, and subsequent 
relaxation led back to the principal (now global) mini- 
mum. 

To gain more insight into the structural relaxations, 
we have found it useful to introduce a measure of the 
"strength of the symmetry breaking." While it is tempt- 
ing to choose a measure that is related to the energy 
difference between the two local minima, like the £ e g in- 
troduced earlier, such a definition has the disadvantage 
of being ill-defined in the single-minimum case. Thus, 
we have adopted instead the following measure. For any 
curve such as that of Figs, l(d-e), 2, or 3, we locate the 
point of minimum d 2 E/du 2 (i.e., vanishing d 3 E/du 3 ) 1 
and then define to be dE/du evaluated at that point. 
We refer to F s b as the "symmetry-breaking force" since it 
measures the strength of the symmetry breaking and has 
units of force. It turns out to have a similar behavior as 
the well-depth difference in the double-well case, but has 
the advantage of remaining well-defined in the single- well 
case. 

Calculations of F s \> in the A-iso and B-iso cases indi- 
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cate that F s b has only a modest and smooth dependence 
on cell volume even when passing through the transi- 
tion from single-well to double-well behavior, confirming 
that the "strength of the symmetry breaking" is not the 
variable parameter in those cases. However, returning 
to the virtual-atom B-hetero system, we find that F s b is 
an extremely sensitive function of 5, with numerical fits 
indicating a 5 3 dependence. It is not hard to see that 
-Fsb must be an odd function of 5, but its cubic behavior 
may seem surprising at first sight. However, the vanish- 
ing of the linear term can be deduced from simple sym- 
metry arguments. An intuitive form of the argument is 
to note that the symmetry-breaking perturbation, which 
is of the form (+5, —5, 0) in successive layers, can be 
regarded as a superposition of two perturbations (26/ 3, 
—(5/3, —5/3) and (5/3, —25/3, 5/3) that do not break 
inversion symmetry. Thus, the principle of superposi- 
tion prevents the occurrence of any symmetry-breaking 
response in linear order in 5, and in particular F s b must 
vanish. (A more systematic analysis may be made by 
considering the Cs v symmetry group consisting of prim- 
itive translations along z and mirrors M z . The ferro- 
electric mode vector and the perturbation 5 are found 
to belong to the A 2 and E respresentations, repsectively, 
and thus cannot couple at linear order.) 

Knowing the form of this extraordinarily strong cu- 
bic dependence of F s b on 5, we can now understand the 
pronounced qualitative differences that were observed for 
the cases of isovalent and heterovalent substitution. We 
find that is about the same, ^0.25mHa/a.u., in the 
B-hetero system with 5 = 0.25 as in the A-iso system. In- 
creasing 5 from 0.25 to 0.4 increases F s b by about a factor 
of four, enough to cause the transition to single-well be- 
havior. A further increase of 5 from 0.4 to 1.0 leads to 
a further increase of the symmetry-breaking force F s b by 
a factor of about 15. Thus, in a fully-developed B-hctcro 
system such as Ba(Sc 1 / 3 Ti 1 / 3 Nb 1 / 3 )03, the strength of 
the symmetry breaking is more than an order of magni- 
tude larger than needed to destroy the secondary mini- 
mum completely. It is hardly surprising, then, that we 
observed no secondary minimum in this case! 

The enormous disparity between the behavior in the 
isovalent and heterovalent cases suggests that it may be 
of great interest to find a way of tuning the system con- 
tinuously from one behavior to the other. One can imag- 
ine doing this by regarding 5 not as a variable atomic 
number, but rather as representing a layer-by-layer com- 
position variable. For example, one could conceive of 
the epitaxial growth of Ba(Sci_ J) Nb l/ )03 in an alternat- 
ing sequence of layers with y = 0.5(1 + 5), y — 0.5, and 
y = 0.5(1 — 5). By the same symmetry arguments, the 
effective strength of the symmetry breaking must again 
scale as 5 3 . Thus, by controlling the concentration vari- 
able 5, one can hope to tune the system over a very wide 
range of behavior. 

One would clearly like to use this tunability to optimize 
the desired characteristics of the material, such as the 



piezoelectric response. It might naively be expected that 
increasing the strength of the inversion symmetry break- 
ing will increase the piezoelectric response. (For example, 
if the parent material at <5=0 were paraelectric, then the 
piezoelectric response would be expected to scale as 5 3 .) 
On the other hand, the materials that have the largest 
piezoelectric coefficients are typically ferroelectrics, and 
since a very strong symmetry breaking suppresses the 
ferroelectric behavior, it might be counterproductive to 
make 5 too large. Clearly, further theoretical investiga- 
tion is needed in order to clarify these issues. 

Many other questions remain open and need to be re- 
solved. For example, the conditions under which the po- 
larization will remain oriented parallel to the growth di- 
rection, the behavior of the system as a function of tem- 
perature, and the properties of materials with simultane- 
ous A-site and B-site substitution, are obvious candidates 
for further study. The thermodynamic behavior of these 
materials, which are technically pyroelectric but interpo- 
late to a ferroelectric limit as 5 — > 0, are also deserving 
of investigation. In the meantime, we hope that our the- 
oretical investigations will stimulate attempts at experi- 
mental growth and characterization of novel perovskites 
with compositionally broken inversion symmetry 
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